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We investigate two closely related subjects: (i) existence of a pure persistent spin current in 
semiconducting mesoscopic device with a spin-orbit interaction's (SOI), and (ii) the definition of 
the spin current in the presence of SOI. Through physical argument from four physical pictures in 
different aspects, we provide strong evidences that the persistent spin current does exist in a device 
with SOI in the absence of any magnetic materials. This persistent spin current is an analog of 
the persistent charge current in a mesoscopic ring threaded by a magnetic flux, and it describes 
the real spin motion and can be measured experimentally. We then investigate the definition of 
the spin current. We point out that (i) the non-zero spin current in the equilibrium SOI device is 
I , the persistent spin current, (ii) the spin current is in general not conserved, and (iii) the Onsager 

relation is violated for the spin transport no matter what definition of the spin current is used. 
These issues, the non-zero spin current in the equilibrium case, the non-conserved spin current, 
and the violation of the Onsager relation, are intrinsic properties of spin transport. We note that 
the conventional definition of the spin current has very clear physical intuition and describes the 
spin motion very well. Therefore we feel that the conventional definition of the spin current makes 
physical sense and there is no need to modify it. In addition, the relationship between the persistent 
spin current and transport spin current, the persistent linear and angular spin currents in the SOI 
region of the hybrid ring, are discussed. Finally, we show that if the spin-spin interaction is included 
into the Hamiltonian, the persistent spin current is automatically conserved using the conventional 
O . definition. 
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I. INTRODUCTION 



In traditional charge-based electronics, the spin degree of freedom has not been fully explored. Until recently, the 
scientists found that spin control and manipulation in the nano-scale can enhance operational speed and integration 
density of conventional charge-based electronic devices. 2,3 In order to successfully use the spin degree of freedom of 
electrons in the conventional semiconductor devices, one has to have a good control, manipulation, and detection 
of the spin and its flow in nano-devices. This emerging field called spintronics is a new sub-discipline of condensed 
matter physics, and it is growing rapidly and generating great interests in recent years. 2,3 

It is well known that the spin-orbit interaction (SOI) plays an important role in the emerging field of semiconductor 
spintronics. SOI couples the spin degree of freedom of electrons to their orbital motions, thereby giving rise to a 
useful way to manipulate and control the spin of electron by an external electric field or a gate voltage. SOI is an 
intrinsic interaction having its origin from relativistic effects that can be quite significant in some semiconductors. 4 
For instance experiments show that the significant SOI indeed exist in some semiconductors, 5-9 e.g. Rashba SOI was 
found in the InGaAs/InAlAs heterostructure or InAs Quantum Wells, 5 ' 6 and both the Rashba and Dresselhaus SOI 
were detected in various III-V material based 2D structures at room temperature, 7 to just name a few. Moreover, the 
strengths of these SOI have been well modulated by the gate voltage experimentally. 8 

Many interesting effects resulting from SOI have been predicted. For example, using the effect of spin precessions 
due to the Rashba SOI, Datta and Das proposed a spin-transistor more than ten years ago. 10 By using the Rashba 
SOI, a method to efficiently control and manipulate the spin of the electron in the quantum dot has been proposed. 11 
Very recently, a very interesting effect, the intrinsic spin Hall effect, is theoretically predicted by Murakami et.al. and 
Sinova et.al. in a Luttinger SOI 3D p-doped semiconductor and a Rashba SOI 2D electron gas, 12,13 respectively, that 
a substantial amount of dissipationless spin current can be generated from the interplay between the electric field and 
the SOI. Since then, the spin Hall effect has generated tremendous interests with a great amount of works focusing in 
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the field of spintronics. 14-19 On experimental side, two groups by Kato et al. 20 and Wunderlich et al. 21 have observed 
the transverse opposite spin accumulations near two edges of their devices when the longitudinal voltage bias is added. 
In additional, a third group by Valenzuela and Tinkham took the electric measurement of the reciprocal spin Hall 
effect, 22 and they have observed an induced transverse voltage in a diffusive metallic conductor when a longitudinal 
net spin current flows through it. 

In this paper, we study two closely related subjects that: (i) we predict another interesting effect that a persistent 
spin current without accompanying charge current exists in a coherent mesoscopic semiconductor ring with symplectic 
symmetry, i.e., with SOI but maintaining the time reversal symmetry; (ii) we examine the issue whether it is needed 
to redefine the conventional spin current, or in other words, whether the conventional definition of the spin current, 
Is = Re{^^vs^} 23 is reasonable in the presence of SOI. Since these two subjects are closely related, we investigate 
them together here. In fact, a few years ago, by using the conventional definition of the spin current Rashba has found 
a nonzero spin current Is in an infinite two-dimensional system with Rashba SOI in the thermodynamic equilibrium. 2 
In his opinion, this spin current Is is not associated with real spin transport, and therefore should be eliminated in 
calculating the transport current by modifying the conventional definition of the spin current. After his work, many 
subsequent work have discussed the definition of the spin current. 25-28 So to uncover the physical meaning of the 
non-zero spin current in the equilibrium system, we must face the question of whether one needs to redefine the spin 
current. 

More than two decades ago, the persistent (charge) current in a mesoscopic ring threaded by a magnetic flux 
has been predicted theoretically, 29 and later observed experimentally in early 1990s. 30 It is now well known that 
the persistent charge current is a pure quantum effect and can sustain without dissipation in the equilibrium case. 
There has also been many investigations on the persistent spin current. 31-35 For example, in a mesoscopic ring with a 
crown-shape inhomogeneous magnetic field 31 or threaded by a magnetic flux 32 , the persistent spin current has been 
predicted and is related to the Berry's phase. Recently, the persistent spin current carried by Bosonic excitations 
has also been predicted in a Heisenberg ring with the magnetic field or in the ferromagnetic material. 33 The reason 
that the persistent spin current exists may be explained as follows. Due to the magnetic field or the magnetic flux, 
there are persistent flows of both spin up and down electrons. In the absence of SOI, this gives rise to the well known 
persistent charge current. In the presence of SOI or magnetic field, the persistent charge current is spin polarized 
resulting a nonzero persistent spin current. Hence the origin of this persistent spin current is the same as that of 
persistent charge current so that the persistent spin current always accompanies with a persistent charge current. 
Recently, 36 ' 37 we have reported that a persistent spin current without accompanying charge current (a pure persistent 
spin current) can sustains in a mesoscopic semiconducting ring with SOI in the absence of the external magnetic 
field, magnetic flux, and the magnetic material. This pure persistent spin current is induced solely by SOI, which is 
different with the pervious case. In main part of this paper, we will give detailed discussions on the existence of the 
persistent spin current and the behavior of the spin current. 

Another motivation of the present paper is to study the definition of the spin current in the presence of SOI. In 
fact, the spin current is the most important physical quantity in the field of the spintronics. So it is important to 
give an appropriate definition of the spin current. This problem, the discussion of the definition of the spin current, 
is first tackled by Rashba. In a recent work by Rashba, 24 he found a nonzero spin current Is in an infinite two- 
dimensional system with Rashba SOI even in the equilibrium case by using the conventional definition of the spin 
current, 23 Is = ReWvsty. So he questioned the correctness of the conventional definition of the spin current and 
suggested that the conventional definition should be modified to eliminate the equilibrium non-zero spin current. 
Besides the nonzero spin current in the equilibrium, there arc others problems with the conventional definition: (1) 
the spin current is not conserved using the conventional definition in the presence of SOI or/and the magnetic field. 
This is because the operator s does not commute with the Hamiltonian H when the system has SOI or/and the 
magnetic field. (2). the Onsager relation is violated using the conventional definition. After Rashba's work, many 
subsequent papers have discussed the definition of the spin current. 25-28 For instance, Sun and Xie suggested that 
in addition to the conventional (linear) spin current I5 = Re^/'vs^f , one needs to introduce the angular spin current 
Io, = Re^ {ds / dt)^ = Re&u) x sVE' to describe the rotational motion (precession) of a spin, 25 because the spin 
operator is a vector. Similar to the conventional linear spin current, the angular spin current can also induce an 
electric field. Shi et. al. 27 gave a new definition of the spin current with I5 = Re^'d(rs)/dt^, in which the operator 
of the spin current Is = d(rs)/dt is the derivative of the whole rs. Compared with the conventional definition of the 
spin current, it has an extra term r(ds/dt) in this new definition. Wang et. al. 28 pointed out that the spin current is 
automatically conserved using the conventional definition of the spin current if the spin-spin interaction is included. 
In the Sec. VI in the present paper, we will explore the definition of the spin current. 

In this paper, we first give physical argument as to why the persistent spin current can exist in equilibrium. In 
order to show that the persistent spin current should exist in the mesoscopic system with solely SOI, four physical 
pictures or physical argument from the different aspects will be discussed: (i) from the picture of electric and magnetic 
correspondence to analyze the driving force of persistent spin current (in Fig. la and Fig. lb); (ii) from the point of view 
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of the spin Berry phase; (iii) from the comparison among four effects, the Hall effect, spin Hall effect, the persistent 
(charge) current, and the persistent spin current (see Fig. 2); (iv) from the point of view of the motion of the spin. 
As an example, we then consider a semiconducting SOI-normal hybrid mesoscopic ring and show that indeed a pure 
persistent spin current can emerge with solely SOI. We note that currently there is no consensus on the definition for 
the spin current in the presence of SOI. In order to avoid the problem of the definition of the spin current, here we first 
use the following approach: we consider a mesoscopic hybrid ring that consists of a Rashba SOI's region and a normal 
region without SOI as shown in Fig.lc. Since there is no spin- flip in this normal region, the definition of the spin 
current in that region is without controversy. So we can calculate and study the spin current in the normal region, and 
to make sure that the persistent spin current indeed can be induced solely by SOI. After making sure the existence 
of the persistent spin current, we then to investigate the definition of the spin current, (i). We point out that the 
non-zero spin current in the equilibrium system in the presence of the SOI is the persistent spin current. It describes 
a real spin motion and has the physical meaning, so this spin current should be kept as it is. (ii). Since spin operator 
is a vector, it has the rotational motion (precession) due to the SOI in accompany with the particle translational 
motion, so in general spin current is not a conserved quantity. In fact, there are experimental indications that the 
spin current is not conserved. 22 On the other hand, in certain cases if one includes the strong spin-spin interaction 
into the Hamiltonian, the spin current obtained from the conventional definition will be conserved automatically, (iii). 
It is well known that the Onsager relation holds under certain condition. We will show that for the spin system, this 
condition does not satisfy. We think that the above three points are actually the intrinsic properties of spin transport 
so we feel that the conventional definition of the spin current need not to be modified. 1 

In addition, we also address the following issues. (1). The relation between the persistent spin current and transport 
spin current. We note that they can not be distinguished from each other in the coherent part of the device. (2). 
In calculating the persistent linear and angular spin currents in the SOI's region of the SOI-normal hybrid ring, 
we find that the persistent spin current still exists in the SOI's region even if the whole ring has the SOI. So the 
normal part is not necessary for generating the persistent spin current. (3). The measurement of the persistent spin 
current is discussed, we suggest that this persistent spin current can be observed by detecting its induced electric 
fields. (4). Including the spin-spin interaction in the ring, we find that the persistent spin current calculated using 
the conventional definition is a conserved quantity. 

The rest of the paper is organized as follows. In Section II, we provide physical arguments and physical pictures 
from four different aspects to show the existence of the persistent spin current. In Section III, we consider a SOI- 
normal hybrid mesoscopic semiconducting ring device to show that indeed a pure persistent spin current can emerge 
in the normal region where the definition of the spin current is without controversy. Then the effect of sharp interface 
between the SOI's region and the normal region, and the relation of the persistent spin current and the transport spin 
current, are discussed in Section IV and V, respectively. In Section VI, we study the definition of the spin current. 
In Section VII, the persistent linear and angular spin current in the SOI's region as well some conserved quantities 
in the hybrid ring device are explored. In Section VIII, we investigate electric fields induced by the persistent spin 
current. This provides a way to detect the persistent spin current. In Section IX, we discuss the effect of spin-spin 
interaction in the ring which makes the persistent spin current conserved. Finally, Section X summarizes the results 
of our work. 



II. PHYSICAL ARGUMENTS FOR THE EXISTENCE OF THE PERSISTENT SPIN CURRENT 

In this section, we argue that the pure persistent spin current (without accompanying charge current) should exist 
in the mesoscopic semiconducting ring device with the SOI. In particularly, this pure persistent spin current can be 
induced solely by SOI even at zero magnetic flux or magnetic field. We examine this new effect from following four 
different aspects. 

A. Analysis of the driving force 

For the persistent charge current in the mesoscopic ring, the magnetic flux or magnetic field acts like a "driving 
force", so one naturally looks for the analogous "driving force" in the spin case. To discuss this question, let us 
consider two devices. The first device consists of a mesoscopic ring (without SOI) where a magnetic atom with a 
magnetic dipole moment is placed at the center of the ring (see Fig. la). In the second device the magnetic atom 
is replaced by a charged atom, e.g., an ion (see Fig. lb). The magnetic atom produces a vector potential A on the 
perimeter of the ring which drives the persistent charge current. By analogy, a charged atom which produces a scalar 
potential </> on the perimeter of the same ring should drive a persistent spin current. 38 Since the presence of this ionic 
center generates a SOI in the relativistic limit, we expect that this SOI which plays the role of the spin "driving force" 
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will induce a pure persistent spin current. In addition, from the picture of electric and magnetic correspondence, the 
persistent charge current in the ring should also change into the persistent spin current, when the center magnetic 
atom (in Fig. la) is substituted by a charged atom (in Fig. lb). 

B. From the point of view of the spin Berry phase 

Second, the existence of the pure persistent spin current can be examined from another point of view using the spin 
Berry phase. 39,40 This physical arguments has been provided in detail in our previous Letter (see the 4-th paragraph 
in Ref. 36 ), so we omit the discussion here. 

C. Comparison among four effects: the Hall effect, spin Hall effect, the persistent (charge) current, and the 

persistent spin current 

In the following, let us compare the four effects: the Hall effect, the spin Hall effect, the persistent charge current, 
and the persistent spin current, from which one expect that the persistent spin current should exist in the mesoscopic 
ring with solely SOI. 

(i) We consider a four-terminal device with a bias V added between the terminals 1 and 3 and a perpendicular 
magnetic field B (as shown in Fig. 2a). For this system there exists the Hall effect, and a charge current is induced 
in the transverse terminals 2 and 4. (ii) Replacing the perpendicular magnetic field by the SOI (see Fig. 2b), a spin 
current emerges (instead of the charge current) in the transversal terminals 2 and 4. This is the spin Hall effect 
predicted recently and generated tremendous interests. 12-18 (iii) Consider a mesoscopic ring with a perpendicular 
magnetic field (see Fig. 2c), a persistent charge current is induced in the ring. This is the persistent current, which 
is well known now. 29,30 (iv) By analogy, a persistent spin current instead of the persistent charge current should be 
induced when the SOI replaces the perpendicular magnetic field (see Fig. 2d). 

Let us discuss the Hamiltonian in four devices in Fig. 2. (i) In Fig. 2a, the Hamiltonian is H = ( - p+ ^^ c - ) — h V(r), 
and there exists the Hall effect because of the vector potential A. (ii) To replace A by the Rashba SOI with 

H = + V(r), the Hall effect change into the spin Hall effect, and a spin current instead of the charge 

current emerges in the terminals 2 and 4 (see Fig. 2b). (iii) In Fig. 2c, the Hamiltonian is H = ( - p+ ^^ c - > — h V(r), and 
there exist the persistent charge current in the ring because of the vector potential A. (iv) By analogy, when A is 
replaced by the Rashba SOI, a persistent spin current should emerge instead of the persistent charge current. 

D. From point of view of the motion of the spin 

By analyzing the motion of the spin, including the translational and rotational motion (precession) of a spin, one 
can also show the existence of the persistent spin current in the mesoscopic device with the SOI. 37 In order to analyze 
the motion of the spin, we need to solve the wave- functions of the electron, so we put this analysis in the appendix, 
where the wave-functions have been solved. 



III. PERSISTENT SPIN CURRENT IN THE NORMAL REGION OF THE SOI-NORMAL HYBRID 

RING 

In this section, we present an example to show that indeed a pure persistent spin current can exist for a mesoscopic 
semiconducting ring with SOI. In the presence of SOI, the spin of an electron experiences a torque and hence Oj 
(i = x, y, z) is not a good quantum number anymore. Because of this, the spin current is not conserved using the 
conventional definition. At present there are controversies on whether one should define a conserved spin current or 
whether there exists a conserved spin current. 25-28 In another word, so far there is no consensus on the definition 
for the spin current in the presence of SOI. In this section, we use the following approach. We first discuss the 
persistent spin current for a one-dimensional mesoscopic semiconducting ring that consists of a Rashba SOI's region 
and a normal region without SOI as shown in Fig.lc. Since there is no spin-flip in the normal region, the spin current 
can be calculated using conventional definition without controversy, so that we can make sure the existence of the 
persistent spin current in the equilibrium case with solely SOI. After we make sure the existence of the persistent spin 
current, we then go back to examine the definition of the spin current, which is given in Section VI. 
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The Hamiltonian of our system is given by: 32,41 



d 2 ia r 



d d 



where E a = ft 2 /2ma 2 , a is the radius of the ring, m is the effective mass of the electron, oy = <j x cos ip + <j y sin</?, and 
n ,p = —a x sin ip + a y cos tp. otfi(ip) is the strength of the Rashba SOI, an{<p) = while < ip < $0j i-e. in the normal 
region, and an(<p) is a constant a R in the SOI's region with $ < <P < 27r. 

The eigenstates of Hamiltonian (1) can be solved numerically in the following way. First in the Rashba SOI's region 
(atR ^ 0), the equation H^(tp) = E^(ip) has four independent solutions %>f°(ip) (i = 1,2,3,4): 32 

*i/2lW - {- S m(9/2)e^+ 1 ^ ) ' (2j 

and ^3/4 = TtyfQ with T being the time-reversal operator. In Eq.(2), the wave vectors k\j 2 = —1/2 + 1/(2 cos 0) ± 
(1/2)^/(1/ cos 2 9) — 1 + 4E/E a , and the angle 9 is given by tan(#) = ctR/(aE a ). Similarly, in the normal region 
(0 < tp < $ ), the Schodinger equation has four independent solutions: *f (ip) = (l,0)^e ikv , ^ 2 (ip) = (1, 0)^ - *^, 
and = T*f /2 , with fc = y/E/E a . Secondly, the eigen wave function ^(<p) with the eigen energy E can be 

represented as: 

di^f(ip), while < ip < $ 
* ( ^ ): > V/ )( >|, >";-.. while * <<p<2ir (3) 

where dj and bi (i = 1,2,3,4) are constants to be determined by the boundary conditions at the interfaces <p — and 
$0- Here the boundary conditions are the continuity of the wave function ^(p)\ ip=0 + = "^(¥>)l v= 27r-/*~ an< ^ the 
continuity of its flux 42 v v ^\ v=Q+ ^+ = v v ^\ v=27r - , where t)^ = (2aE a /ih){d/dip + (i/2)a r tan(#)] is the velocity 
operator. By using the boundary conditions, we obtain eight series of linear equations: 

ai + a 2 + cos(0/2)e i27rfel &i + cos(6/2)e i27rk2 b 2 

+ sin(6»/2)e- i27r ' £1 6 3 + sin(6i/2)e- j27rfe2 64 = (4) 
a 3 + a 4 - sin(0/2)e i27rfel &i - sin((9/2)e j27rfe2 6 2 

+ cos(6»/2)e~ j27r ' £l 6 3 + cos(6>/2)e- j27rfe2 64 = (5) 
e ife *°ai + e^ fe *°a 2 + cos(0/2)e ifcl *°&i + cos^Je** 2 * ^ 

+ sin(0/2)e- J ( fel+1 )* o 63 + sm(9/2)e'^ k2+1 ^% 4 = (6) 
e~ ik9 °a 3 + e ik ^ a 4 - sm(9/2)e l< - kl+1 ^°b 1 - sm(6>/2)e i ( fe2+1 > < H 2 

+ cos(9/2)e- tkl ' s '% 3 + cos{9/2)e- tk2 ' i ' a b i = (7) 
fcoi - fca 2 + ^(Aije* 2 '* 1 ^ + A x (k 2 )e l2 ^b 2 

-^ 2 (A:i)e- j2rfl 63 - A 2 (fc 2 )e- j2T H 4 = (8) 
-ka 3 + ka 4 - A 2 (k 1 )e i2wkl b 1 - A 2 {k 2 )e l2 ^b 2 

-A 1 {k 1 )e~ i2 ^b 3 - A 1 (k 2 )e- t27rk2 b 4 = (9) 
fce ife *°ai - ke- lk * a a 2 + A 1 {k 1 )e lk ^ a b 1 + A x (k 2 )e lk ^ a b 2 

-i4 2 (fci)e- i ( fcl+1 )* 6 3 - A 2 (/c 2 )e- J ( fe2+1 ) <I ' 6 4 = (10) 
-ke- ik *°a 3 + fce lfe *°a 4 - A^e^+^h - A 2 (k 2 )e l ^ +1 ^°b 2 

-Ai(fci)e- ifcl *°&3 - ^(fc^e"** ^ = (11) 
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where A x (x) = x cos(0/2) - \ tan(0) sin(0/2) and A 2 (a;) = (a;+ 1) sin(0/2) - \ tan(0) cos(0/2). The eigenvalue E n can 
be solved numerically by setting the determinant of the coefficient of the variables and in the above eight scries 
of linear equations to zero. 

Now we present the numerical results. Fig. 3a shows the eigen values E n versus the Rashba SOI's strength a_R. For 
the normal ring (an = 0), the eigenvalues arc n 2 E a with fourfold degeneracy, and the corresponding eigenstates are 
(1, 0)^e ±m ' f and (0, i)t e ±»«y. 4 3 As the SOI is turned on the degenerate energy levels split while maintaining twofold 
Kramers degeneracy. The higher the energy level, the larger this energy split. Typically, the splits are on the order 
of E a at an = 10 _11 eV r m, with E a w QA2meV for the ring's radius a = 50nm and the effective mass m = 0.036m e . 
The eigenvalues E n versus the normal region's angle $o are also shown (see Fig. 3b). For $o = 27r, the whole ring is 
normal and E n are fourfold degenerate. 43 When <J> is away from 27r, the degenerated levels are split into two, and 
the splits are larger with the smaller $ . When $ = 0, the whole ring has the Rashba SOI, and the split reaches the 
maximum. 

Since E n is twofold degenerate, we obtain two eigenstates for each E n , which are labeled V& n (^) and T^P n ((p) . 44 
With the wave functions, the spin current contributed from the level n can be calculated straightforwardly using the 
conventional definition Ig^ip) = Re^^v^a^ n {i — x,y,z). Notice that the spin current Isji is a tensor, where 
i, j = r,<f,z in the cylinder coordinates and i, j — x,y,z in the orthogonal coordinates. The first index j describes the 
direction of the motion of the electron, and the second index i represents the direct of the spin. Because the device 
in the present paper is a ring, the motion of the electron has no components along radial (r) and z axis. Hence the 
spin currents Ig ri and Ig zi ( Is( r /z)i^P) = ^- e ^n^r/z^i^n) are zero. Only spin current Ig vi with the electron moving 
along the if direction is non-zero. To simplify the notation, hereafter we use the symbol 72^ to replace Ig vi - Since 
there is a controversy about the definition of spin current in the SOI's region, we will calculate the spin current only 
in the normal region in this section. The spin current in the SOI's region will be studied in the section VII after the 
definition of the spin current is investigated in the section VI. In the normal region, the spin current is conserved, so 
Isi(<P) 1S independent of the angle coordinate (p. 

Fig. 4 shows the spin current Ig i versus the Rashba SOI's strength ur for $ = 7r. 45 Since E n is twofold degenerate, 
the wave- functions can be arbitrary combination of ty n (ip) and T^ n ((/?). But the spin current remains the same. 
Our results in Fig. 4 show that Ig x is exactly zero for all level n while $o = and Ig and Ig z exhibit the oscillatory 



pattern with a^. A 7r/2-phase shift between Ig y and Ig z is observed with . / {Ig y ) 2 + (7g z ) 2 approximately constant. 



For two adjacent levels 2n — 1 and 2n, their spin current have opposite signs, and 7|" _1 + 7|" = if ap> = 0. We 
note that the spin current Ig t is quite large. For example, the value E a is equivalent to the spin current of a moving 
electron in the ring with its speed 4 x 10 5 m/s. 

The spin current 7^ versus the angle $o that describes the normal region at a fixed ajj = 3x 10~ 11 eVm is shown in 
Fig. 5. When $o = 27r, the whole ring is normal, Ig y and Ig z are exactly zero. However, Ig x is non-zero (see Fig. 5a). 
Note that 7 2 ,™ -1 +7|™ = (n = 1, 2, 3, ...) and I Sx = at <&o = 27r, so the total spin current 75^ is still identically zero 
because that the (2n-l)-th and the 2n-th states are degenerate and have the same occupied probability at $o = 27r. 
When <!>o ^ 2tt where part of the ring has the SOI, three components Ig x / y / z of the spin current can be non-zero. 
For the larger n, the absolute value of the spin current |7^J is larger. For two adjacent levels 2n — 1 and 2n, their 
spin current have opposite signs, which is similar to the result of Fig. 4. 

Now we calculate the equilibrium total spin current I$i contributed from all occupied energy levels: Isi — 
2 Y,„ I Stf( E n), where f(E) = l/{exp[(E - E F )/k B T] + 1} is the Fermi distribution with the Fermi energy E F 
and the temperature T and the factor 2 is due to the Kramers degeneracy. The persistent charge current and the 
equilibrium spin accumulation are found to be zero because the system has the time-reversal symmetry. Fig. 6a, b 
show the total spin currents Isi versus the Rashba SOI's strength an for different Fermi energy Ep. One of the main 
results is that the spin current indeed is non-zero when ^ 0. The persistent spin currents Igi in Fig. 6 have the 
following features. At aji ~ the whole ring is normal, so Isi is exactly zero. With increasing I$i increases 
initially and then oscillates for the large ckr. In Fig. 6a, b, the parameter $o is 7r, i.e. half of the ring is normal and the 
other half of the ring has the SOI, then Is x is zero, and only I$ y and 7g z are non-zero. If $o 7r, the components 
Isx/y/z can be non-zero. At certain there is a jump in the curve of Isi versus a^. This is because for this 
the Fermi energy Ep is in line with a level E n , leading to a change of its occupation. At zero temperature, the jump 
is abrupt as shown in Fig.6a,b. But at finite temperature, this jump will be washed out. In fact, these results are 
similar to the persistent (charge) current in the mesoscopic ring. 29 

The spin current Isi versus the angle $ of normal region at a fixed = 3 x IQ^^eVm is shown in Fig. 6c, d. 
When <!>o — 27r, the whole ring is normal, we have Isx/y/z — 0- When <I>o is away from 27r, the spin current Isx/y/z 
emerges. For some Fermi energy a jump appears in the curve Isi-^o (as shown in Fig. 6c), which behaviors is similar 
as the jump in the curve Isi-ctR- For other Fermi energies, however, the jump in the curve Isi versus $o (see Fig.6d) 
disappears when the Fermi energy Ep is not in line with the level E n at all values of $0- In particular, in the limit 
when $o goes to zero, i.e. when there is no normal region in the ring, the spin current 75^ and Ig z still exist. This 
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indicates that the normal region is not necessary for generating Isi- 

In the above numerical calculation, the temperature T is set to zero. Now we consider the effect of the finite 
temperature T. Fig. 7 shows the persistent spin currents Is y /z versus the strength of SOI a R at $o = t with the 
different temperatures ksT. When the temperature ksT = 0, the spin currents \Isi\ are the largest and 1st shows a 
jump if the Fermi energy Ep is in line with a level E n . With the increase of the temperature from zero, this jump 
is smoothed, and the spin currents \Isi\ at these a R near the point of the jump (e.g. 7.4 x 10 _11 eUm in Fig. 7) 
decreases sharply even at the very low temperature. But the spin currents \Isi\ at these a R that are far away from 
this point of jump (e.g. a R < 5 x 10 _11 eUm in Fig. 7) are not effected very much by the temperature fc^T even for 
ksT reaching 0.2E a . Upon further the temperature ksT is raised and on the same order of the energy-level interval, 
the spin currents 1st is reduced for all values a R . This is because the probability of occupation, i.e. f(E n ), of the 
energy level vary smoothly versus the level index n and the spin current due to the adjacent levels are opposite in 
sign. However, even when fcgT reaching ksT = 0.5E a , the persistent spin current Isi is still quite large. For the 
ring's radius a = 50nm and the effective mass m = 0.036m e , E a is about 0A2meV. Then the temperature T is 
approximately 2.1K at ksT — 0.5E a . This temperature T is easily reached at the present technology. 8,46 



IV. DISCUSSION THE EFFECT OF SHARP INTERFACE BETWEEN THE SOI REGION AND THE 

NORMAL REGION 



In the above section, the coefficient Rashba SOI a R ((p) vary sharply in the interface of the normal and SOI's 
part. Now we examine the effect of sharp interface. Let us consider a hybrid ring device with the SOI coefficient 
a R (<p) varying continuously along the ring. The Hamiltonian is same as Eq.(l), with a R (p) = for < ip < tt and 
a R (<p) = a R sin 2 (tp) for n < tp < 2n. In this case, both a R (ip) and da R {tp) / dp are continuous at the interfaces ip = 
and tt. In this section, we show that persistent spin current still exists for a hybrid ring device with SOI varying 
continuously along the ring. 

For this continuous varying Rashba SOI coefficient a R (tp), the Hamiltonian can not be solved analyti- 
cally. Here we numerically solve this Hamiltonian by using the discrete tight-binding model. Notice that 



2a 



discrete becomes, 47 



ia R (<p)a v /(2a) = -± a R (<p)a r -g- + -g-a r a R (p) 



and the discretized Hamiltonian 



AT 



JV 



f * a rJ a R,j + °"r,j + l a fl,j + l 



^ +1 2a 



2Aip 



dj + H.c. 



(12) 



where cij and are annihilation and creation operators at the point j, 



N is the number of the points in the ring, 



'r,j 



Vr(j&V) = &x COs(jAp) + 



a v sm(jAip), 



Ap — 2tt/N is the angle between two neighboring points, t — E a j (Ap) 2 

and a R .j — aji(jAcp). In the above Hamiltonian, the point index iV + 1 is same with the point index 1. Then by 
calculating the eigen- values and eigen- vectors of the Hamiltonian matrix with the dimension N, the eigen- values E n 
and the eigen wave-functions ^f n j (f nj = ^ n (jApj) of the ring device can be easily solved. 44 After solving the eigen 

wave-functions f nj , the spin current Ig k (ip) — Re^\v v a\^ n — j^R e {— U^n^fe^*" = x iVt z ) m the normal 
region can be obtained straightforwardly from 



l Sk,j 



2Aip 



(13) 



In order to show that the above method is correct and reliable, we first solve the model of the section III again, with 
the sharp varying SOI's coefficient a R (<p): a R (ip) — while < ip < $0 and a R (ip) — a R while $ < ¥ < 2?r. The 
results of the eigen-energies are shown in Fig. 8. When the number of points N = 20, we can see that the eigen-energies 
from the above discrete method are quite different from the exact values obtained from the method in the section III 
(see Fig. 8a). But with increasing TV (e.g. N = 50), this difference become very small (see Fig. 8b)). When TV = 150, 
the eigen-energies from the discrete method are in excellent agreement with the exact results (see Fig. 8c). This means 
that our results using the above discrete method converges for large N. 

Since E n is twofold degenerate, the arbitrary combination ci$> n {<p>) + c 2 T 1 $ l n ((p) still is the eigen wave- functions, 
so the wave-function is un-certain. In the following we examine the correctness of our spin current Ig i (i — x,y, z), 
which is calculated from the wave- functions. Fig. 9 shows the spin currents Ig contributed from the level n for the 
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case of the sharp varying SOI, and it shows that the results of the spin current of the above discrete method are also 
in excellent agreement with the exact result at large N. In particular, it is surprising that the number of points iV 
need not to be very large. For N — 50, the difference between numerical and the exact results is already very small, 
and for N = 150, there is almost no difference. 

Now we are ready to examine the effect of sharp interface. Fig. 10 shows the results of the hybrid ring with SOI 
varying continuously along the ring, with an{ip) = for < <p < n and otn(ip) = cxr sin 2 (<p) for it < tp < 2tt otherwise. 
Fig. 10a, b, and c are for the eigen values, the spin current Jg contributed from the level n, and the equilibrium total 
spin current I$i in the normal region versus cur, respectively. The eigen values are four-fold degenerate at cur = 0, 
and this degeneracies are split into two twofold Kramers degenerated states. The spin current Ig contributed from 
the level n oscillates with a R = 0. These results of the eigen values and Ig are similar to case of the sharp interface. 
In particularly, the results show that the (equilibrium) persistent spin current Igi is still non-zero and has a quite 
large value (see Fig. 10c, d). This indicates that the persistent spin current I$i indeed is originated from the SOI, and 
it is not the artifact of sharp interface. 

V. RELATION BETWEEN THE PERSISTENT SPIN CURRENT AND THE TRANSPORT SPIN 

CURRENT 

Through the physical arguments and physical pictures from four different aspects in the section II, and the analytic 
results of an example of a SOI-normal ring (in which the definition of the spin current is without controversy in the 
normal region) in the section III, as well as the discussion concerning the interface in the section IV, so far we have 
plenty evidence to demonstrate that the persistent spin current indeed exists in a mesoscopic semiconducting ring 
device with an intrinsic SOI. This persistent spin current can be induced solely by a SOI, and it exists in an equilibrium 
mesoscopic device without a magnetic field, a magnetic flux, and in the absence of any magnetic materials. Besides 
the above mentioned ring geometry, we find that the device can also be of other shapes, 37 ' 48 e.g., a disc device, a 
quantum wire, a two-dimensional system, etc. Thus, it is a generic feature that a pure persistent spin current appears 
in a system with SOI. In this section, we will discuss the relation between the persistent spin current and the normal 
transport spin current. 

First let us review and discuss the relation of the persistent (charge) current and the transport (charge) current 
in a mesoscopic ring device. To consider a mesoscopic ring threaded by a magnetic flux coupled to two leads (the 
lcad-L and the lead-R) which act as the electron reservoirs (see Fig. 11a). There exists dissipation in the two leads 
and they are always in equilibrium for an isolated lead-L(R). The size of the ring is assumed within the coherent 
length and no dissipation in the ring. At zero bias (i.e. the equilibrium case), the transport charge current is zero 
everywhere, including the two leads and the ring, but the persistent charge current exists only in the ring, but not 
the leads because of presence of dissipation in the leads. If a non-zero bias is added between two leads, a transport 
charge current flows from one lead through the ring to the other lead, so there exists both transport and persistent 
charge currents in the ring but only the transport current in the leads because of dissipation. In other words, at a 
finite bias, the charge currents iei(e2) m the two arms of the ring are the sum of the transport currents I e i(e2),t and 
the persistent current I e ^ p : I e \ — I e \,t + Ie,p and I e 2 — I e 2.t — Ie,p- Can one distinguish the transport charge currents 
Iei(e2).t and the persistent charge current J eiP ? In fact, they (l e u e 2),t an d I e ,p) can not be distinguished either in 
theory or in experiment. The transport charge currents I e i( e 2),t and the persistent charge current I eiP in the ring 
have identical behaviors, both of them are dissipationless 49 , capable of inducing a magnetic field, etc. So in principle 
only the total charge currents / e i(e2) m the arms are observable physical quantities. If some dissipative impurities 
are introduced in the ring, then the transport charge currents I e u e 2),t show dissipation while the persistent charge 
current I e , p does not. Under this circumstance, can one distinguish the transport currents / e i(e2),t and the persistent 
currents I e . p ? It turns out that one can still not distinguish these types of currents. When dissipative impurities 
arc introduced, the original dissipationless ring is changed into a new different ring, in which the persistent charge 
current I ep is quenched while the transport charge currents I e i( e 2),t are normally reduced due to dissipation. Thus, 
it is impossible to obtain I e i(e2),t and I e;p of the original dissipationless ring. 

The relation between the persistent spin current and the transport spin current is identical to the relation between 
the two charge currents discussed above. Consider that part of a mesoscopic ring contains a SOI but without a 
magnetic flux, and with two leads coupled to this ring (see Fig. lib). In the equilibrium case, a persistent spin current 
emerges in the ring, neither the transport spin current nor the persistent spin current is present in the leads. Under a 
spin-motive force, 19 ' 50,51 a transport spin current flows from one lead through the ring to the other lead. The persistent 
spin current in the leads is always absent independent of with or without a spin-motive force due to dissipation. On 
the other hand, with a spin motive force, both persistent spin current and the transport spin current exist in the ring. 
Similarly to the charge currents, both persistent spin current I SjP and transport spin currents I s i( S 2),t in the two arms 
of the ring are indistinguishable since they behave identically in all physically measurable properties. Both of them 
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are dissipationless 49 , describing the real spin motion, capable of inducing an electric field, and so on. 

VI. THE DEFINITION OF THE SPIN CURRENT 

From the sections II, III, and IV, we have made sure that the persistent spin current exists in the equilibrium 
mesoscopic semiconducting device in the presence of SOI. In this section, we examine the definition of the spin current. 
The first work to question the conventional definition of the spin current is by Rashba. 24 After that, many subsequent 
papers have discussed the definition of the spin current as mentioned in the introduction. 25-28 In summary, 27 in the 
presence of SOI one faces three problems when using the conventional definition Ig — Re^vs^-?^ (i) there exists a 
non-zero spin current even in the equilibrium system, (ii) the spin current usually is non-conservative, and (iii) the 
Onsager relation is violated. Therefore suggestions have been made in previous papers that one needs to modify the 
conventional definition of the spin current. In the following, we examine these three problems in detail and we argue 
that there is no need to modify this conventional definition 23 Ig = Re^vs^l. 

(i) the non-zero spin current in the equilibrium system: 

From the investigation and the discussion in the sections II, III, and IV, we have clearly shown that this non-zero 
spin current is the persistent spin current. It describes the real motion of spins, has the physical meaning, and can 
be observed in the experiment in principle (see the section VIII). So this non-zero spin current should be kept in the 
calculation of the total spin current. 

(ii) non-conservation of the spin current: 

In this sub-section, we argue that in general the spin current is not conserved in the presence of the SOI or/and 
the magnetic field. However, in certain cases, the spin current can be conserved by including the strong spin-spin 
interaction. Let us analyze this problem from the both aspects of the theory and experiment. In the aspect of the 
theory, we give two examples of the non-conserved spin current (for intuition the readers can also consider the spin 
classic vector). 

First, we want to show that in general if the system has spin flip mechanism, the spin current will not be conserved. 
For instance, in the presence of a rotating magnetic field or circular polarized light 52 the spin current exists in a 
device connected to only one terminal. This one terminal device and its spin translational motion and precession are 
as shown in Fig. 12a. Here at x < 0, the spin is without precession, and the spin pointing to the y-direction moves 
along the +x-direction while the spin pointing to the — y-direction moves along the opposite direction. So it has a 
non-zero spin current I atXy at x < 0. Near the point x = where the quantum dot located, the spin precesses and 
is flipped due to the presence of a rotating magnetic field or circularly polarized light. 52 Then the spin accumulation 
does not vary with time and the system maintains the steady state. It is obvious that the spin current is not conserved 
because of only one terminal in this device. 

In the second theoretical example, to consider the spin translational motion and precession on a ring as shown in 
Fig. 12b. At point A with its angular coordinate ip = 0, the spin pointing to +x-dircction moves down while the spin 
pointing to — x-direction moves up, then the non-zero element of the spin current is Is,yx- At another point B with 
its angular coordinate ip = n/2, the spin pointing to +j/-direction moves along the +x-axis and the spin pointing to 
— y-direction moves along the — x-axis, then the non-zero element of the spin current is Is.xy So it is obvious that 
the spin current is not conserved in this ring device, but the spin accumulation still keep invariant. In fact, this 
example is similar to the persistent spin current in the mesoscopic ring while the whole ring has the SOI (see the 
section VII. C and the appendix), in which the motion of the x-y- plane elements of the spin is as in the Fig. 12b and 
the spin accumulation is zero everywhere. 

In addition, there arc experimental indications that the spin current is not conserved. For example, let us consider 
the experimental result by Valenzuela and Tinkham. 22 In this experiment, they have clearly shown that a pure spin 
current injects from the FM1 electrode into the Al strip, reduces with its flowing forward because of the spin flip, and 
finally disappears while the distance much longer than the spin diffusion length (see Fig. la, b, c, and Fig. 4 in Ref. 22 ). 
This experimental result gives a strong proof of that the spin current should be non-conserved in the presence of the 
SOI, the magnetic field, the magnetic impurity, or others. 

(iii) violation of the Onsager relation: 

The Onsager reciprocal relation is an important theorem of the near-equilibrium transport theory. Up to now, 
the Onsager reciprocal relation is always satisfied for transport of any physical quantity (e.g. the charge transport, 
thermal transport, etc) by suitably defining a corresponding current. However, in the following we point out that the 
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Onsager relation for the spin transport is in general violated. In particular, it is impossible to restore the Onsager 
relation regardless of how to modify the definition of the spin current. This is very different from all previous cases. 
In this sub-section, we first recall the tenable condition of the Onsager reciprocal relation, and then examine the case 
of the spin current. We find that the spin transport does not always meet this condition no matter how to define 
the spin current, although this condition is met for all previous transport phenomena people have studied. So the 
Onsager relation is in general violated for the spin transport. 

First, let us recall the tenable condition of the Onsager relation. Considering the currents / = {Ij} and its 
corresponding forces F = and they have the following relationships: 



ijFj, (14) 

where is the conductivity. If the local entropy production dS/dt per unit time can be expressed as: dS/dt = 
IiFi, there exists an Onsager relation Gij = Gji (assuming that the system has the time-reversal symmetry). 
Next, we give an intuitive example to show the tenable condition of the Onsager relation. Let us consider the charge 
conductivity in the two-dimensional system. To take the vectors e x and e y (as shown in Fig. 13a) as the base vectors 
of the charge current and its force [i.e. the gradient of the potential, W(r)], the equation (14) becomes: 



h \ = ( G xx G xy \ ( V X V 
ly J \ Gyx G yy J \ VyV 



(15) 



In this case, the Onsager relation is tenable, and G xy = G yx = G . However, if taking the non-orthonormal vectors c x 
and e r (as shown in Fig. 13b) as the base vectors of the current and the force, the Onsager relation will be violated as 
shown in the following. For the base vectors e x and e r , any current vector I and force vector W can still be expressed 
as: I — I x e x + I r e r and VV = V x Ve x + \7 r Ve r . Then it is easily to obtain the relation between (V X V, ^ y V) and 
(V x V,V r V): 

(J'V)=V(Z'YA, and(l*Y r \=^ (W * V 

where 



VyV J \ V r V J ' V \7 r V J \ V y V 



sin tp 

So under the base vectors e x and e r , the Eq.(15) changes to: 
Ix \ tt-1 I G xx G \ TT / \7 X V 



u=(irr) as) 



i, J u [gT (;"„ :i ) 11 { v,i 

G xx - G cos ip I sin tp (G xx - G yy ) cos tp - G a cos(2</?)/ sin tp \ / V X V 
G / sin ip Gyy + G cos ip/ simp J I V r V 



(17) 



It is obvious that the two off-diagonal elements of the conductivity in the above equation are not equal, so the Onsager 
relation is violated in the non-orthonormal base vectors e x and e r . In fact, for the charge current, the Onsager relation 
is only tenable under the orthonormal and linear independent base vectors. 

Now, let us discuss the spin current. The spin current has 3x3 = 9 elements and the charge current has 3 elements. 
So here the current / (including spin and charge) and the corresponding force F totally has 12 elements and the 
conductivity has 12 x 12 = 144 elements. For simplicity and clarity, we consider a one-dimensional system and the 
electron can only move along the x axis. In this case, the current / and the force F have only 4 non-zero elements, 
and they are: 

I (ls XX 7 Is X y : Is X Z : Ic X ) ) (1^) 

F = {V x V sx ,V x V sy ,V x V sz ,V x V c ), (19) 

where V s i (i — x, y, z) is the spin chemical potential and V c is the (charge) chemical potential. 19,50,51 The spin chemical 
potential V s i means that the electron of the spin pointing to +i and — i-direction occupy up to V c + V S i/2 and V c — V s i/2 
(shown in Fig. 13c, d, e), respectively. Then the Eq.(14) becomes: 



( ^sxx N 




( G xx 


G xy 


G xz 


Isxy 


- 


G yx 


Gyy 


Gyy. 


^sxz 




G zx 


Gzy 


G zz 


\Icx J 




V G cx 


G C y 


G cz 



G X c 

Gyc 

G zc 



( v x v sx 

V ' x^sy 

V X V SZ 



(20) 



G cc / V V X V^ 
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and the conductivity has 4 x 4 = 16 elements. 

Now we show that the spin transport does not meet the tenable condition of the Onsager relation. First, the three 
base vectors Oi (i — x, y, z) are not orthonormal, and it is impossible to find a series of the orthonormal base vectors 
regardless of what combination of the three <7j. 

In addition, once the spin chemical potential of one component is fixed, the other two spin potentials V S i can not 
exist. For example, if we give the value for the spin potential V sz , this means that the electron of the spin along 
+z direction occupies up to V sz /2 and the electron of the spin at — z direction occupies up to —V sz /2 if V c = (see 
Fig.l3e), 19 ' 50 and the electron occupational state has completely been determined. Then one can not further specify 
the spin states and the corresponding occupation number along +x and —x (or +y and —y). Therefore, while giving 
the value for one component of the spin potential, e.g. V sz , the others V sx and V sy do not exist. This conclusion hold 
even for spin free and conserved system with [a t , H] — for all i (i — x, y, z). Due to the fact that three spin potentials 
V S i can not be simultaneously evaluated, the Eq.(20) [i.e. the Eq.(14)] does not exist for the spin transport, regardless 
of how to define the spin current. Therefore there is no Onsager relation in spin transport, and it is impossible to 
restore the Onsager relation by modifying the definition of the spin current. In fact, the Onsager relation is not 
satisfied for any existing definitions of the spin current. 

Furthermore, we can also use the four normal terminal device (as the device in Ref. 19 ) to examine the Onsager 
relation. Because the four terminal leads are normal metal without the SOI, there is no spin-flip. For this system, 
there is no controversy for the definition of the spin current in the terminal leads. Therefore we can avoid the definition 
of the spin current and examine the Onsager relation. The results also show that the Onsager relation does not exist 
for the spin transport. This gives an additional proof. 

We wish to mention that because V sz and V c (or V ax and V c , or V sy and V c ) can be determined simultaneously, the 
Onsager relation for G zc — G cz (or G xc = G cx , or G yc = G cy ) might exist, 19,27 e.g., in the suitable boundary condition 
in the four terminal device. 19 However, it is impossible that six relations (G zc = G cz , G xc — G cx , G yc — G cy , G xy = 
G yxi G xz = G zx , and G yz — G zy ) are satisfied simultaneously. 19 For three-dimensional systems, the conductivity 
has 12 x 12 = 144 elements Gij, similar conclusion applies, i.e., it is impossible that Gij = Gji are satisfied for all 
off-diagonal matrix element simultaneously. 

(iv) discussions 

From the discussion in the last three sub-sections, we have clearly shown that (i) the non-zero spin current in the 
equilibrium SOI's device is the persistent spin current, (ii) in general the spin current is not conserved, and (iii) the 
Onsager relation is violated for the spin transport. In particular, it can not be restored regardless of how to modify 
the definition of the spin current. Therefore the three "flaws" of the conventional definition of the spin current, which 
has been mentioned and commented in some previous papers, are intrinsic properties of spin transport. In addition, 
the conventional definition has very clear physical intuition and has described the spin motion very well. Using this 
conventional definition, one can account for many effects that relate the spin transport, 25 e.g., the heat produced 
by the spin current, the spin currents induced electric field, and the force and the torque acting on the spin current 
in the presence of electric field. Therefore we make the conclusion that there is no need to modify the conventional 
definition I s = Re^^vs^. 23 

Before the end of this section, we wish to mention that if the spin-spin interaction is included into the Hamiltonian, 
the spin current calculated using the conventional definition will be conserved automatically. 28 We will discuss this 
in detail in section IX. We also note that our discussion above does not contradict with the angular spin current in 
our previous paper by Sun and Xie. 25 To see this, let us first recall the (linear) velocity v and the angular velocity 
uj. (a) The (linear) velocity v = dr/dt (or the velocity operator v = dr/dt) can describe the translational motion of 
the vector (or the rigid body) and there is no need to modify this velocity definition v = dr/dt. (b) On the other 
hand, the vector has the rotational degrees of freedom except for its translational motion, so we need to introduce the 
angular velocity Co to describe its rotational motion. We emphasize that two statements (a) and (b) do not contradict 
to each other and can be rephrased as (a) The (linear) spin current I s = Re&vsW can describe the translational 
motion of the spin and there is no need to modify its definition, (b) On the other hand, the spin has the rotational 
degrees of freedom (precession) except for its translational motion, so we need to introduce the angular spin current 
I w = Re^(ds/dt)^f = Re^^Co x to describe its rotational motion. 53 

VII. THE PERSISTENT SPIN CURRENT IN THE REGION WITH SOI 

After clarifying the definition of the spin current, we return to discuss the persistent spin current in the SOI's region 
of the normal-SOI ring in this section. In the sub-section A, we investigate the linear and angular persistent spin 
current. In the sub-section B, we show that the persistent spin current still exists even using the new definition of the 
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spin current as in Ref.[27]. The case of entire ring with SOI is studied in the sub-section C. At the last sub-section 
D, we present some conserved quantities that are related to the persistent spin current. 



A. the linear and angular persistent spin currents in the SOI's region 



From the definition of the (linear) spin current Igi — Re^ '(1/2) (v v Si + Sji^)*, for a ring device we have: 



' Isx(<fi) = Im{^[a x aE a £ 



dip 



^cos^]*} 
^Msin^]*} 



Isyfa) = Im{^[a v aE a -L , 2 
[ I St (<p) = Im{&[* t aE a £]*} 

To make the above equations discrete, Isi(<p) changes into Isij — Isi(jAtp): 



(21) 



L Sy,j 



l Sz,j 



= aE a Im < * 
= aE„Im\^ 



2E, 



2E, 



i sin(jAp)tf j 



aE, 



(22) 



Let us calculate the angular spin current I w , 25 which describes the rotational motion (precession) of the spin. For a 
ring device with the Hamiltonian of Eq.(l), one has: 



dsx/dt = (— ia z j 2a) {a R(tp) sin tp,d/ dip}, 
d§y/dt = (ia z /2a){an(ip) cos(p,d/dip}, 
ds z /dt = (-i/2a){a R (tp)a (p ,d/d(p}, 



(23) 



where {A, B} = AB + BA. Using the cylindrical coordinates, ds/dt — (ds r /dt,ds v /dt,ds z /dt) = 
^{Q,—(j z {aR{ip)-,d/dip},{aii{ip)a tp ,d/dip}). Then the angular spin current can be calculated from its definition, 25 

I w = Re^(ds/dt)^ straightforwardly, so does its discrete version: 



2Aip 



+ 



2Aip 



(24) 



Since the wave function ^ n {f) for each eigenstate is known from the section III or IV, the (linear) spin current Ig i 
and the angular spin current 7" 4 from the n-th state can be easily obtained. Using the exact method in the section 
III or the discrete method in the section IV, we find that the results for and I" i match perfectly. Fig. 14 shows the 
linear spin current Ig x i y / Z (i.e. the element in the orthogonal coordinates) and the angular spin current I™ v / Z (i.e. the 
element in the cylindrical coordinates) versus the angular coordinates tp. In the normal region, Ig x i y / Z is conserved 
and is independent of the coordinates <p, and I™ v / Z is zero, since the spin has only the translational motion without 
the precession there. On the other hand, in the SOI region, except for J" r = 0, the (linear) spin current Ig X / y / z ana - 
the angular spin current I™ v i z are all non-zero and non-conserved, because of the spin precession in the presence of 
the SOI. The linear spin current Isx/y/ziv) versus <P is always continuous, even at the interface between the normal 
region and the SOI region where the strength of SOI azn(ip) changes abruptly. But the angular spin current 7™ v ^(<p) 
versus <p> shows a jump whenever an(ip) versus ip has as abrupt change, as shown in Fig.l4d, e, and f. The jump 
position of 7™ , is located at the abrupt point of a R (ip). The spin currents Jg f (and 7") versus ip for the states with 
the same-parity, e.g. n = 1, 3, 5, or n = 0, 2, 4, 6, are similar in shape (see Fig. 14a, c, d, and f), but the value of 
l-Tsil ( or i s mucn larger for a larger n. Thus, the spin current from the highest occupied level dominates in the 
(total) persistent spin current. 

Fig. 15 shows the total persistent spin current Isx/y/z versus the angular coordinates <p for different Fermi energies 
Ep. The persistent spin current is conserved in the normal region but not so in the SOI region because of the spin 
precession. Due to the fact that the spin current l^'fy^ from the highest occupied level maxn dominates in Isx/y/zi 
Isx/y/z behaves similarly as 7g^™ z (see Fig. 14 and Fig. 15). 
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B. calculating the spin current using the definition in Ref. [27] 

We have clarified and demonstrated that the conventional definition of the spin current makes sense in the section VI. 
However, in this sub-section, we present our calculated results of the persistent spin current by using the new definition 
of the spin current as appeared in Rcf.27, and discuss its consequences. In the normal region, the persistent spin current 
Is is the same regardless which definition is used. But in the SOI region, I5 depends on the definition. In the following 
let us discuss I5 in the SOI region. By using the new definition in Ref. 27 , I 5 = Re^ = R e ^ [{,§ + vds/dt] 
the spin current from the level n, then the total persistent spin current, depends on the choice of the coordinate origin. 
First, if setting the origin at the center of the ring, the element Isipi (i = x,y,z) with the spin motion along the ip 
direction is Re^^v^s^ . So Is V i is completely the same to the result using the conventional definition, and it is still 
non-zero in the equilibrium case and non-conversed in the presence of a SOI. On the other hand, the element Isri 
(i = x,y, z) with the spin motion along the radial direction is aRe^ (dSi/ dt)^ = al^i, so it is also non-zero, but the 
same element is exactly zero using the conventional definition. The element Ig zi {i — x, y, z) with the spin motion 
along the z direction is zero, same as with the conventional definition. Second, if the coordinate origin is not located 
at the center of the ring, all 9 elements of the spin current are in general different to those by using the conventional 
definition. In particular, they are all non-zero and non-conserved in equilibrium. 

C. the case when the entire ring is with a SOI 

Let us consider the case of $0 — > 0, i.e. the normal region is gradually getting smaller and at the end the whole 
ring contains the SOI. Fig. 16 shows the persistent spin current Isx/y/z versus the angle coordinate ip for $0 = tt/2, 
7r/4, and 0, respectively. The results clearly show that the persistent spin current Isx/y/z does exist, and its value 
\Isx/ y /z\ is even larger with decreasing of the normal region, i.e. <j> . Eventually when the entire ring contains the 
SOI, \Isx/y/z\ reaches its maximum value. This means that the normal region is not necessary for the existence of 
the persistent spin current. 

In fact, if the whole ring has a constant SOI with an(ip) = a^, the persistent spin current can be analytically 
obtained. In this case, the eigenwave function is: 32 

where n = 0, ±1, ±2, ... and the eigenvalue E n is given by 

E n = E a [n 2 + (n + 1/2) (1 - l/cos0)] (26) 

Then the linear spin current lg and the angular spin current I™ from the state n are as follows: 

Ig = —E a F{6) [e x sin # cos + e y sin sin 93 - e z cosO] (27) 
I™ = (E a /a)F{9) [e x sin 9 sin <p — e y sin 9 cos <p] 

= - (E a / a) F(9)e v sin 9 (28) 

where F{9) = [n + 1/2 — 1/(2 cos 6)]/2tt. The persistent linear and angular spin currents I5 and \ u are obtained by 
summing Ig and I™ over the occupied states. In addition, from the wave function [Eq.(25)] and the spin currents 
[Eq. (27,28)], the spin motion in the ring can be obtained (see the discussion in Appendix). 

D. three conserved quantities in a ring device 

From the results of Eq. (27,28), three quantities characterizing the spin current lg are found to be conserved, 
although the spin current itself is not conserved in the presence of a SOI. (a) The spin current l Sz with spin po- 
larization along z-direction is a conserved quantity for the ring geometry, (b) The magnitude of the spin current 

Tg = ^J{Isx) 2 + i^Sy) 2 + (^Sz) 2 = E a F(9) is a constant of motion, (c) For a given cross-section of the ring, the vector 
of spin polarization for the spin current Tg makes a fixed angle with the normal direction of that cross-section. This 
angle is a constant for any cross-section of the ring. In this sense, the spin current Ig is "conserved", although the 
direction of spin polarization for Ig is not a constant of motion due to a SOI. So the non-conservation of spin current 
in the ring device means that while moving along the ring the direction of spin polarization is precessing due to the 
torque from the SOI. For the hybrid ring, our numerical results also show that the magnitude of spin current Ig is 
again a constant of motion across the hybrid ring, but the element of the z-direction Tg z and the angle in (c) are not. 
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VIII. THE INDUCED ELECTRIC FIELD BY A PERSISTENT SPIN CURRENT 

There are a number of experiments that have been carried out to confirm the existence of spin current, 20-22 e.g. to 
observe the spin current induced spin accumulations by the Kerr effect, 20 or to make the electric measurement through 
the reciprocal spin Hall effect. 22 Since the persistent spin current is an equilibrium property, the above mentioned 
methods are not suitable. There is also a proposal that a spin current may cause a spin torque that can be measured 
experimentally 54 ' 55 . Very recently, Sonin pointed out that this method can be employed to detect the persistent spin 
current. 55,56 On the other hand, we note that the persistent charge current can be detected by measuring its induced 
magnetic held. 30 ft has been shown that the persistent spin current can also generate an electric field. 25 ' 33,57 ' 58 So 
this offers another way to detect the persistent spin current by measuring its induced electric field. In the following, 
we calculate the persistent spin current induced electric field and electric potential, and show that this electric field 
or the electric potential can be observed in the present technology. 

The induced electric fields Es( r ) an d E w {v) at space point r = (x,y, z) by the linear and angular spin currents I5 
and I w in the ring device are: 25 



where [1b is the Bohr magneton and r' (ip) = a(cos (p, sin ip, 0) is the position vector in the ring. Considering the whole 
ring having a constant SOI, the persistent linear and angular spin currents have been solved in the section VII. C [see 
the Eqs. (27,28)]. Substituting them into the above formulas, Eqs. (29,30) (note that Isi is the element Is v % of the 
linear spin current) and with the help of the first and second kind complete elliptic integral functions K(x) and E(x): 

K{x) = / -x(smtp) 2 dip, (31) 

Jo 

E(x) = / VI - x(smcp) 2 dtp, (32) 
Jo 

the induced electric fields Es( r ) and E u (r) can be obtained straightforwardly. The electric fields Es( r ) an d E u (t) 
are rotational invariant about the z axis, and in the plane e z -e r , i.e. the elements Es<p and E uv are zero. So here we 
only show Es( r ) and E u (r) in the x-z plane with r = (x, y, z) = (x, 0, z): 

E Sx = 2 3 C {K(A)R\ [a (R 2 + R 2 _ - z 2 (a 2 + x 2 + z 2 )) cos9-z (R 2 + R 2 _ + a 2 (a 2 - x 2 + z 2 )) sintf] 

CLXix ix j 

—E(A) [a (R 2 + R 2 _(a 2 + x 2 ) - 16a 2 x 2 z 2 ) cos (9 - z (R 2 + R 2 _(2a 2 + x 2 + z 2 ) + 8a 2 x 2 (a 2 - x 2 - z 2 )) sin6»] } (33) 
E Sz = f 4 {K(A)R% [az(-a 2 + x 2 + z 2 ) cosO + (R 2 + R 2 _ + 2a 2 z 2 ) sinO] 

CLi\ Jl j 

+E(A) [az (-R 2 + R 2 _ + 8a 2 (a 2 - x 2 + z 2 )) cosO + (R 2 + R 2 _(3a 2 - x 2 - z 2 ) - 8a 2 z 2 (a 2 + x 2 + z 2 )) sin0] } (34) 

E - = I a^Rj [( ° 2 +X2 + Z2)E{A) R + K{A)] (35) 
K. - [(« 2 - - 2 - * 2 )E{A) + R\K{A)] (36) 

and Es y = E uy = 0. Here R\ = (a±x) 2 + z 2 , A = —Aax/R 2 _, and c = —^,Qg^BE a F{6)/h. Then the total electric field 
Et = Eg + E W are also easily obtained, and Et can be expressed as a gradient of a potential V(r), Er(r) = — VV(r), 
where 

V(r) = 2C 2 [R 2 + K(A)cos9+ ((a 2 - x 2 - z 2 ) cos - 2az sin 9) E(A)] . (37) 
R—R_l_ 



In fact, this total electric field Et can also be expressed as: 

(,n\ . 

\r-v'(tp)\ 



E T = -W = -cV / P e (p) • „,; J^ adip, (38) 
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i.e. is equivalent to the electric field generated by a ID electric dipolc moment P e (f) = 

— (cos 9 cost/?, cos 9 sin <p, sin 9) = —e v x (— sin#e r + cos9e z ) oc v x S in the ring (see Fig. 17c). 

Fig. 17 shows the electric-field lines of E s , E w , and E T in the x-z plane. The across points of the ring and the 
x-z plane are at (a, 0,0) and (—a, 0,0). The electric-field lines have the characteristics: The field lines are in the x-z 
plane and E$ y — E^ y — Et v = 0. V x E5 and V x E w usually are non-zero, but the total electric fields Ex has the 
behavior V x Et = 0, i.e. § Et • dl = 0. The electric fields in the x-z plane are mirror symmetry around z axis with 
Es/w/t,x(x,z) = -E s/u / T , x (-x,z) and E S / u / T ,z( x > z ) = e s/u>/t,z(- x > z )- 

Fig. 18 shows the electric-field strengths E Sx / z , E ux / Z , and E Tx / z , along the two horizontal dashed lines [from the 
point (0,0, a) to the point (4a, 0, a), or from the point (0,0,0.1a) to the point (2a, 0,0.1a)] in the Fig. 17. In this 
calculation, we consider that only the lowest level n = (the ground state) in the ring device is occupied by the 
electron, i.e. taking the parameter E < E F < E Xl and in this case I5 = l° s and I w = 1° . We also take the parameters 
g factor g = 2 and the efficient electron mass m = 0.036m e . At the point x = 0, the x-direction element Es/ u /t,x = 
due to the mirror symmetry, but the ^-direction element Es/ u /t,z keep quite large value still. In slightly far away 
from the ring device (e.g. z — a in Fig. 18a, c, e), the electric field E w induced from the persistent angular spin current 
is in the same order with the electric field Eg induced from the persistent linear spin current. So E w is important 
for contributing to the total electric field E<p. On the other hand, while very near the ring device (e.g. z = 0.1a in 
Fig. 18b, d, f), Es is much larger than E w and E5 is dominant in Et. It is worth to mention that the total electric 
field Et can reach 10~ 2 V/m at the point (a, 0, 0.1a) which is 0.1a = bnm over the ring (a, 0, 0) (see Fig.l8f). Also let 
us estimate the electric potential difference due to Et, this potential difference between two points (a, 0,0.01a) and 
(a, 0, 0.01a) is about InV. Although this potential value is very small, it is measurable in the present technology. 25 ' 33 



IX. THE SPIN-SPIN INTERACTION AND THE CONSERVED PERSISTENT SPIN CURRENT 

In previous sections single spin picture is adapted, i.e., there is no spin-spin interaction, we find that in general the 
spin current is not conserved. In this section, we demonstrate that if one includes strong spin-spin interaction, the 
spin current will be conserved. In this case, if a spin precesses, the response of other spins will be against it, i.e, they 
will precesses in opposite directions. As a result the total spin precession is zero everywhere and therefore the spin 
current is automatically conserved by using the conventional definition. 28 ' 59 To account for the spin-spin interaction, 
the Hamiltonian H is: 

H = J2 Ho(r t ) + J2 J(r t ,r 3 )S(r t ) . S(r,), (39) 

i i,j 

where Hq(t) is the one-body Hamiltonian (e.g. the Hamiltonian in Eq.(l) for the ring device) and the second term is 
the spin-spin interaction. Usually, it is very difficult to solve this Hamiltonian because of the many-body interaction, 
J{ r i, r .?')S(ri)»S(rj). By introducing an induced self-consistent field Hi that could be due to spin-spin interaction, 28 

the many-body Hamiltonian H in Eq.(39) reduces into the one-body form: 

ff(r) = ff (r)+<7«i?i(r) (40) 

Now the Hamiltonian of Eq. (40) is easily solved, and the linear spin current and the spin torque all depend on H\ . 
The induced self-consistent field Hi in Eq.(40) is determined by requiring that the spin torque (or the angular spin 
current) is zero or l w (Hi) = for any r. Once Hi is solved the persistent spin current is automatically conserved by 
using the conventional definition. 

Now we apply this method to the ring without the normal region and set SOI to a constant with an(<p) = an- We 
start with the following Hamiltonian 

H = H + a r Hi (41) 

where H is the original Hamiltonian [see Eq.(l)] for a ring with full SOI. The reason that we choose the induced field 
as a r Hi is because we know I z is conserved and there is a torque along e v . So a self-consistent induced magnetic 
field is needed (a term a r Hi) to balance the torque. The eigenfunction of the new Hamiltonian is the same as Eq.(25) 
given by 

T , , f cost 9/2) e mv \ 
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where the eigenvalue looks the same as before 

E n = E a [n 2 + (n+ 1/2) (1 - I / cos 6)] (43) 

but 9 takes a different value, 

tanfl = -^r + - - Hl — (44) 
aE a E a (n+ 1/2) K 1 

Note that if the self-consistent field Hi is zero, tan# recovers the non-interacting case. Using the conventional 
definition the persistent spin current with polarization in three directions are easily calculated and found to be 

I Sx = -E a Fi{6) sin 9 cos tp 
Is y = — E a F\{9) s\n9svcnp 

hz = — E a [ti + 1/2 — 1/(2 cos 9)\ cos 9 (45) 

where 

Fi{6) = [n + 1/2 - a/(2aE a sm6)]/2TT (46) 

If Hi = 0, Fi(9) = [n + 1/2 — 1/(2 cos 9)]/2n = F(9) that is the previous result without the spin-spin torque 
interaction. Now we calculate the total torque. From continuity equation, the torque is just the angular spin current 
l w = Re^u x sVE'] with uj = H\e r — (ia/a)d/dip e v x e z . It is easy to show that I w has only e v component 

I uv =Hi cos 9+-(n cos 9 -sin 2 9/2) (47) 

Note that the purpose of introducing the self-consistent field Hi is to make sure that the spin current is conserved or 
the torque vanishes as we have discussed earlier. Setting the torque I utp to zero we obtain the second equation that 
determines the self-consistent field Hi 

Hi = ^— (n cos 9 -sin 2 9/2) 

a cos 9 

= ~{n +1/2 -1/(2 cos 0)) (48) 

Once Hi is solved from Eqs.(44) and (48), the spin current will be conserved following the continuity equation. 
Plugging Eq.(48) into Eq.(44), we find 

Oi 

(n+1/2 sin0= — — (49) 
laE a 

This means that Fi{9) defined in Eq.(46) is zero and the conserved persistent spin current is nonzero only for the spin 
polarization along z-direction 

Isz = -E a [n + 1/2- 1/(2 cos 9)] cos 9 (50) 

where 9 is determined by Eq.(49). 

We wish to point out that even for the hybrid ring, the above approach can be used and the persistent spin current 
is also conserved in the presence of spin-spin interaction. The only difference is that we have to introduce three 
self-consistent fields Hi, i = 1,2,3. From the energy dispersion relation, we have relationship between 9 and Hi. By 
requiring the torque along each direction to be zero, we obtain three additional equations. These four equations will 
determine Hi and 9. This in turn gives the displacement spin current and hence the conserved persistent spin current. 



X. CONCLUSION 



In summary, we have investigated two closely related subjects: (a) the prediction of a pure persistent spin current in 
an equilibrium mesoscopic device with solely spin-orbit interaction (SOI), and (b) the issues concerning the definition 
of the spin current. Through the physical arguments and physical pictures from four different aspects, the analytic 
calculation results of a SOI-normal hybrid ring, as well as the discussion of the sharp interface between the normal 
and SOI part, we demonstrated that the persistent spin current indeed exists in the equilibrium device with a SOI 
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alone. In particularly, we emphasize that this persistent spin current is an analog of the persistent charge current 
in the mesoscopic ring threaded by a magnetic flux, and it describes the real spin motion and are experimentally 
measurable. 

After showing the existence of the persistent spin current, we investigate the definition of the spin current. We 
point that: (i) the non-zero spin current in the equilibrium SOI's device is the persistent spin current; (ii) in general 
the spin current is not conserved; and (iii) the Onsager relation is violated for the spin transport, and in particularly, 
it can not be recovered through modification of the definition of the spin current. So these three "flaws" , the non-zero 
spin current in the equilibrium case, the non-conserved spin current, and the violation of the Onsager relation, of the 
conventional definition of the spin current are intrinsic properties of spin transport. In particular, the conventional 
definition, 1$ = Re^'vs^, possesses a very clear physical picture, and is capable of describing the spin motion. So we 
draw the conclusion that the conventional definition of the spin current makes physical sense, and no need to modify 
it. 

In addition, a number of problems have also been discussed. The relation between the persistent spin current and 
transport spin current is discussed, and we find that they are indistinguishable in the coherent part of the device. 
We calculate the persistent linear and angular spin current in the SOI's region of the hybrid ring, and the results 
show that the persistent spin current still exists in the SOI's region, even when the SOI covers the whole ring. The 
measurement issue of the persistent spin current is also discussed, we suggest that the persistent spin current can be 
observed by detecting its induced electric field. In the presence of a spin-spin interaction in the ring, we find that the 
persistent spin current using the conventional definition is automatically conserved. 
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Appendix 

In this appendix, we analyze the motion of the spin which relates to the persistent spin current in the equilibrium. 
For simplicity, we consider the constant SOI case, and the eigen wave- function has been solved in Eq.(25). From 
this wave-function, the spin S = (\I/„|s|\E'„) = |(— sin 9 cos ip, — sin 9 sin ip, cos 9). This spin vector is in the e z -e r 
plane, and its angle with the e z (i.e. z) axis is —9. Then the spin motion can also be obtained straightforwardly by 
solving the velocity and its angular velocity w. The direction of spin (translational) motion is counter-clockwise, while 
precessing with uj in the perpendicular direction of 5*, so that the spin is in the e z -e r plane all along. The element 
of the spin in the x-y plane and its motion are shown as in Fig. 12b. On the other hand, for the time-reversal state 
T\I>„, the spin direction, its motion direction, and the precession direction, all reversed, as shown the clockwise arrow 
in Fig. 12b. But the spin current of T\&„ is completely same with that of \I/„. Therefore the persistent spin current 
indeed describes the real motion of the spin. 

In fact, besides the ring geometry, the device also can be other shape. For example, we have analyzed the spin 
motion and the persistent spin current in the quasi one-dimensional equilibrium quantum wire. 37 Similar conclusions 
can be drawn. 
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FIG. 1: (Color online) (a) and (b) are the schematic diagrams for a mesoscopic ring with a magnetic atom or an ion at its 
center, (c) Schematic diagram for a hybrid mesoscopic ring having the Rashba SOI in part of the ring while the other part 
being normal. 

FIG. 2: (Color online) (a), (b), (c), and (d) are schematic diagrams for the devices of the Hall effect, the spin Hall effect, the 
persistent (charge) current, and the persistent spin current, respectively. 

FIG. 3: (a) and (b) show the eigen energies E n vs. or for $o = it and vs. <3>o for or = 3 x 1CT 11 eVm, respectively. The radius 
of the ring a — 50nm. 

FIG. 4: Ig y (a) and Ig z (b) vs or for $o = it and a = 50nm. Along the arrow direction, n = 7, 5, 3, 1, 0, 2, 4, 6, and 8. Here 
the level index n = 0, 1, 2, represent the ground state, the first excited state, the second excited state, respectively. 

FIG. 5: I% x (a), (b), and I£ x (c) vs. $ for a R = 3 x lO^eVm and a = 50 urn. Along the arrow direction, n — 6, 4, 2, 0, 
1, 3, and 5. Here the level index n = 0, 1, 2, represent the ground state, the first excited state, the second excited state, 
respectively. 

FIG. 6: (Color online) (a) and (b) show Is y and Isz vs aR for <l?o = it- (c) and (d) show I,sx/ y /z vs the angle of the normal 
region 3>o for oh = 3x 10~ n eVm and Ef = 3E a (c) and 6E a (d). The radius of the ring a = 50nm and temperature T = 0. 

FIG. 7: (Color online) (a) and (b) show I$ y and Is z vs a R at the different temperature T. The radius of the ring a = 50nm 
and $o = it- 

FIG. 8: (Color online) The eigen-energies E n vs. cur. The red dotted curves are for the tight-binding model of the discrete 
Hamiltonian with the lattice points N = 20 (a), N = 50 (b), and N — 150 (c). The black solid curves are from the exact 
method used in the section III, i.e. these black solid are completely the same as the curves in Fig. 3a. All parameters are the 
same as in Fig. 3a. 

FIG. 9: (Color online) I$ y vs. a«. The red dotted curves are for the tight-binding model of the discrete Hamiltonian with the 
lattice points N = 20 (a), N = 50 (b), and N = 150 (c). The black solid curves are for the exact method used in the section 
III, i.e. these black solid are completely the same as the curves in Fig. 4a. All parameters are the same as in Fig.4a. 

FIG. 10: (Color online) (a) the eigen energies E n vs. aR, (b) Ig y vs. or, (c) the persistent spin current Is y vs. or, and (d) 
the persistent spin current Isz vs. aj for the ring with a non-sharp interface (see the text). The ring radius a = 50nm and 
the number of lattice points N = 200. In (b), along the arrow direction, n = 7, 5, 3, 1, 0, 2, 4, 6, and 8. In (c) and (d), the 
temperature T = 0. 

FIG. 11: (Color online) (a) is a schematic diagram for the ring device threaded by a magnetic flux and coupled to left and 
right leads, (b) is a schematic diagram for the ring device with the SOI and coupled to left and right leads. 

FIG. 12: (Color online) (a) Schematic diagram for a spin movement in a terminal device, in which the spin moves along the x 
axis, and flips about x = 0. (b) Schematic diagram for the spin movement and precession in a ring. 

FIG. 13: (Color online) (a) and (b) are schematic diagrams for the base vectors e x and e y , and e x and e r . (c), (d), and (e) are 
schematic diagrams for the spin potential V sx , V sy , and V sz , respectively. 

FIG. 14: (Color online) The linear spin current Is x / y / z and the angular spin current I" v / Z vs the angle coordinates ip for the 
level n — (a and d), 1 (b and e), and 2 (c and f). The parameters are aR — 3 x 10~ n eVm, $o = it, and the ring radius 
a = 50nm. The solid (dark) curve, dashed (blue) curve, and dotted (red) curve in panels (a, b, and c) correspond to Ig x , 
Is y , and Is z , respectively. The solid (dark) curve and dotted (red) curve in panels (d, f, and e) correspond to IZ V and 7" z , 
respectively. 

FIG. 15: (Color online) The persistent spin currents Isx/y/z vs - the angle coordinates ip for the Fermi energy Ef = 2E a (a) 
and 6E a (b). The temperature T — and the other parameters are same with Fig. 14. The solid (dark) curve, dashed-dotted 
(blue) curve, and dotted (red) curve correspond to Is x , Is y , and Is z , respectively. 

FIG. 16: (Color online) The persistent spin currents Is x / y / z vs. the angle coordinates ip for <E>o = ir/2 (a), 7r/4 (b), and (c). 
The temperature T = and the other parameters are same with Fig. 14. The solid (dark) curve, dashed-dotted (blue) curve, 
and dotted (red) curve correspond to Is x , Is y , and Is z , respectively. 

FIG. 17: (Color online) Schematic plots of electric-field lines of the electric fields E u (a), Es (b), and Et (c) for the case when 
the entire ring has a constant or. 

FIG. 18: (Color online) The electric field strengthes E UJX i z (a and b), Eg x i z (c and d), and E Tx / z (e and f) vs. the position x 
along the horizontal dashed lines in Fig. 17 with z = a (a, c, e) and z — 0.1a (b, d, f). The parameters are aR = 3x 10~ n eVm, 
<&o = 0, the ring radius a = 50nm, and only the lowest level n — occupied. The solid curve and the dotted curve correspond 
to E w/S/T „ and E u/S /t,z, respectively. 
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